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Abstract 

Non-adiabatic contribution of environmental degrees of freedom yields effective 
inertia of spin in effective spin dynamics. In this paper, we study several aspects of 
the inertia of spin in metallic ferromagnets. (i) a concrete expression of the spin inertia 
ms'- ms = hSc/{2gsd), where Sc is the spin polarization of conduction electrons and 
Qsd is the sd coupling constant, (ii) dynamical behavior of spin with inertia, discussed 
from viewpoints of a spinning top and of a particle on a sphere, (iii) behavior of spin 
waves and domain walls in the presence of inertia, and behavior of spin with inertia 
in the case of a time-dependent magnetic field. 


1 Introduction 

Different from other fundamental quantities such as mass and charge, spin is a dynamical 
quantity, and its dynamics have been widely studied and applied in science and technology. 
In particular, recent rapid growth of spintronics provides a stage where deeper understand¬ 
ings of spin dynamics directly lead to practical applications. 

The dynamics of spin are governed by the spin Berry phase [1] , and its equation of motion 
includes only the hrst-order time derivative of spin. This is natural because spin is an angular 
momentum and its equation of motion takes the familiar form: the time-derivative of the 
angular momentum (i.e. spin) is given by the torque acting on it. Without any torque, the 
solution of the equation of motion of spin is only a static one. This is in contrast with, for 
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example, the case of a massive point particle, which has inertia and can move at non-zero 
speed as its free motion. In this sense, spin does not have inertia. 

However, for systems where spin interacts with other environmental degrees of freedom, 
spin dynamics are affected by those environmental degrees of freedom, and the dynamical law 
of spin is changed to be an effective one. For example, in metallic ferromagnets, conduction 
electrons affect the dynamics of localized spins (i.e. spins of atoms on lattice cites). A typical 
effect is spin-damping (e.g. 0 ). where the energy and the angular momentum of spins are 
transferred to the environmental degrees of freedom and, as the result, spins relax to their 
ground state within a certain time scale. Without this effect, spins undergo the Larmor 
precession around the applied magnetic held forever. In this respect, the existence of the 
degrees of freedom other than spins changes the dynamical behavior of spins signihcantly. 
In the equation of motion of spins, this damping is represented by the Gilbert damping term 
[3]. This damping term includes only the hrst-order time derivative of spins. Therefore, spin 
still does not have inertia even when we take into account the Gilbert damping effect. 

The effects of the environments other than the Gilbert damping can be studied system¬ 
atically by the derivative expansion, where the effects are expanded in powers of the time 
(and spatial) derivative of spin. From that point of view, the Gilbert damping term gives 
the leading order term in that expansion. In the higher orders, there appear terms which 
include the second-order time derivative, the third-order time derivative and so on, in the 
equation of motion of spin. These terms with higher order time derivatives are interesting 
in that, like the Gilbert damping term, they change the dynamical law of spin itself, more 
than give additional torque on spin. In particular, from a comparison with the form of the 
Newton’s equation of motion of a massive point particle, the term with the second-order 
time derivative of spin plays the role of the inertia of spin. 

Such spin inertia has been discussed in the literature (e.g. ref. |1]). In particular, recent 
progress in ultrafast magnetization mu motivated several works. In refs. [Hi, the inertia 
of spin was introduced phenomenologically and was shown to give additional nutation to the 
motion of spin. The time scale where the effect of the inertia is signihcant was discussed, 
based on the work of Brown [3], to be sub-picosecond order. In ref. HDl, the equivalence 
between the dynamics of spin with inertia and a spinning top was discussed. Microscopic 
derivation of spin inertia was performed in refs. m and [12] . In ref. m, an extended breath¬ 
ing Fermi surface model was used, and the relation between the Gilbert damping coefficient 
and the spin inertia was given in terms of physical quantities (such as Fermi-Dirac occu¬ 
pation numbers) of conduction electrons. The time scale for the nutational motion to be 
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damped by the Gilbert damping was estimated to be sub-picosecond order. In ref. [12], a 
general expression of the contribution of the conduction electrons to spin dynamics was 
discussed. The spin effective dynamics was shown to be non-local in general, which can be 
approximated as local dynamics by the derivative expansion of spin. The inertial term of 
spin arises in that derivative expansion and its general expression was given in term of the 
Green’s function of the conduction electrons. 

Since the inertia of spin is conceptually interesting in its own and gives the hrst step 
toward the understanding of non-adiabatic contribution of environmental degrees of freedom 
to spin effective dynamics, further investigations are worthwhile. Although the expressions 
of the inertia of spin were given as integral forms mm, an explicit expression of the inertia 
of spin in terms of parameters of a model has not been obtained so far. Furthermore, the 
effect of the inertia of spin on the dynamics of spin has been discussed only for spatially 
homogeneous spin system under time-independent magnetic held. In this paper, we present 
detailed theoretical study of the effects induced by the spin inertia based on an sd model. 
In section [H we derive a concrete expression of spin inertia in terms of the parameters in the 
sd model. In section [3l the basic behavior of spin with hnite inertia is studied with the help 
of its two equivalents: a symmetric spinning top and a massive charged particle on a sphere 
subject to a monopole held. In section H] we study spatially inhomogeneous system, and 
discuss that spin waves and magnetic domain walls acquire an additional oscillation mode 
due to spin inertia. We also study the behavior of spin under large and time-dependent 
magnetic held, and hnd an unusual behavior of spin where the velocity of spin is parallel to 
the direction of the time-derivative of magnetic held. 

2 Spin effective action and inertia 

The inertia of spin arises naturally in its ehective dynamics, which takes into account the 
ehects of the environmental degrees of freedom [T2|. Consider a system where a classical held 
of localized spin, S{x,t) = Sin{x,t) (with Si = [S'! hxed), and a held c{x,t) representing 
the other environmental degrees of freedom are interacting with each other. For concrete¬ 
ness, we consider the case of metallic ferromagnets in this paper. Conduction electrons are 
represented by annihilation and creation operators, c and c. The total action is given by 
5s [n] -|- 5e[c, c, n] where 5s [n] is the action of spin n and 5e[c, c, n] is that of conduction 
electrons with their interaction with spin n. When we are interested only in the dynamics 
of spin n, it is convenient to integrate out c and c, and derive the ehective action of spin n. 
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The contribution from electrons is given by path integration as 


exp 



A5efr[n] 


J VcDc exp 


-5e[c, c, n] 


( 2 . 1 ) 


The sum of this AiSeff[n] and the original spin action iSs[n] gives the total spin effective 
action. 

It is difficult to calculate AiSeff[n] exactly, so we should rely on a perturbative analysis. 
We here perform derivative expansion, where AiSefr['/^] is expanded in powers of d^n (/r = 
t,x,y,z). When the system is isotropic, the general form is 


A5eff[n] = 


d^x 


dt 


Sc<j){cos6 — 1) — 


+ 0{{d^nf), 


( 2 . 2 ) 


where n = (sin 6* cos 0, sin 6^ sin 0, cos 6*) and i = x,y,z. We have divided the Lagrangian 
density entirely by lattice volume so that each coefficient represents a quantity per each 
lattice cite. The hrst term is the spin Berry phase with Sc the spin polarization of the con¬ 
duction electrons and the second is the spin-spin exchange interaction induced by electrons 
with Jc the coupling constant. In the hnal term, there arises the inertial term of spin with 
the inertia This rUg has the dimension of [kg ■ m^], the same as that of the moment of 
inertia. In ref. [12], a general expression of spin inertia is derived in term of the Green’s 
function of the conduction electrons. Let us here calculate a concrete expression of spin 
inertia. As in ref. [12], a typical example of iSe[c, c, n] is the sd model [13], where conduction 
electrons interact with localized spins n as 


5e[c, c, n] = / d'^xdt c [ ihdt -|- 


2m 


+ Cf + gsdfi ■ cr 


(2,3) 


Here, m is the mass of the conduction electrons, is the Fermi energy, is the sd coupling 
constant, and cr is the Pauli matrix vector. To obtain the derivative expansion of AiSefr[n], 
we perform SU{2) gauge transformation c —)■ U{x, t)c with an SU{2) matrix U acting on the 
spinor indices, so that the sd interaction becomes diagonal, c(n ■ a) c ^ ca^c. Due to this 
unitary transformation, there appears a so-called spin gauge held = —iU^d^U in (12.3p 
through d^c —)■ U{d^ -|- iA^)c. This contains the hrst-order derivative d^n. Therefore, 
expanding exp(iiSe[c, c, n]//l) in powers of A^, we can calculate AiSeff[n] perturbatively in 
powers of d^n. The spin polarization Sc and the inertia can be calculated as (see 
Appendix |A| for details) 


fc r,3 _ ^3 

Sc = a^- ^^ 

2 67r2 


rUg = 


hSc 

25^sd 


(2.4) 
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with {2m) = ep ± g^d- The inertia can be rewritten as 


ms = {kFa)^-^—f{gsd/eF), f{x) = ^ [(1 + x) 2 - (1 - x) 2 
OTT^ ep 6X L 


(2.5) 


For 0 < X < 1, the function /(x) is a only slightly decreasing function from /(O) = 1 to 
/(I) ■= 0.94. Therefore, spin inertia does not depend much on the sd coupling constant 
gsd and is proportional to when a and kp = yj2mep/h are regarded as independent 

parameters. 

Adding this AiSeff[n] to the original spin action of the form 

d^x 


5s nl = 


-dt 


Si(j){cos9 — 1) — 


( 2 . 6 ) 


with Ji the exchange coupling between spins, we obtain the total spin effective action as 
5efr [n] = Ss [n] + A5efr [n]. Including the Zeeman coupling with external magnetic helcj^ we 
obtain 


Sp.ffln] = 


d^x 


dt 


/■ ^2 ™ 
SB ■ n + S^{cos6 — 1)- —{di'n)'^ + — 


5 . 2 

n 


+ 0{{d,n)% (2.7) 


with S = Sc + Si the total spin amplitude per lattice cite and J = {JcS^ + JiSf)/S‘^. We 
have set the gyromagnetic ratio as unity. 

As we will see below, spin with hnite inertia has a typical precession mode with frequency 
Uq ~ S/mg- Using mg fl2.4p or 02.51) . and assuming S h, kpa 71 and ei? ~ leV, the 
energy scale of this frequency becomes JtjjJq ~ leV, so that its period is 271/ojo ~ O.lps. 
Therefore, as far as this simple estimation suggests, the existence of the inertia is signihcant 
for the dynamics of sub-picosecond scale. 

The equation of motion derived from this effective action 02.7p is 


Sh = —SB X n — JS'^d^n x n + rngii x n. 


( 2 . 8 ) 


Thus, the inertial term produces acceleration-dependent torque. We can rewrite this equa¬ 
tion of motion, by taking vector product with n, as 

rngii = Sn x h + SB + JS'^dfn — {SB ■ n + rngh^ — JS‘^{di'n)‘^)n. (2.9) 

The Gilbert damping effect adds a term —aSh, with a the dimensionless constant, to the 
right hand side of the equation of motion 02.9p . Therefore, the Gilbert damping plays the 

^To be precise, there emerge other coupling terms between spin and electromagnetic field in addition to 
the Zeeman coupling, such as in ref.[T3], by integrating out conduction electrons. In this paper, we simply 
assume that they are negligible. 
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same role as the familiar linear damping force for a point particle, and the time scale for 
this damping term to be signihcant is tdamp ~ ms/{o.S). On the other hand, the time scale 
for the inertial term to be effective is, as we will see below, tinertia ~ rris/S. Therefore, for 
a 1, we can neglect the Gilbert damping term as long as we are interested in the dynamics 
within the time scale tmertia- 

The equation of motion fl2.8p can be rewritten (when S = 0) in the conservation form 
of the angular momentum current (j°, J*), 

+ dij^ = 0, where = Sn + rngn x n, = JS'^diU x n (2-10) 


(5o = d/dt). Note that the angular momentum (per lattice cite), which is the Noether 
charge corresponding to the invariance of the action fl2.7p under S'0(3) rotation in the in¬ 
ternal spin space (see Appendix [B] for details), is no longer proportional to n but includes 
the non-adiabatic contribution of the conduction electrons, m<jn x n. Originally, the total 
angular momentum consists of that of the localized spin and that of the conduction elec¬ 
trons: = Sin + {ha?/2) {ccrc) with Si the amplitude of the localized spin. In the lowest 

order of the derivative expansion, i.e. in the adiabatic limit, the spin of the conduction 
electron aligns with that of the localized spin, so that {ho?/2){ccrc) = Sen with Sc the spin 
polarization of the conduction electrons. Beyond the adiabatic limit, the direction of the 
spin of the conduction electron is generally different from that of the localized spin. The 
derivative expansion incorporates this difference systematically, and the next order term in 
{ho? / 2 ){ c ( tc ) is given by m^n x n. 

Remarkably, the relation between the inertia and the spin polarization Sc in fl2.4p . 
rris = hSc/{2gsd), is easily obtained without any microscopic calculation, as follows. As we 
have discussed in the last paragraph, the angular momentum Aj^ derived from AiSefr[n] 
represents the spin polarization of the conduction electrons. 


= Sen + rusn x n = — — {cere). 


ho? 


Since {ho? /2){c(tc) obeys the following equation of motion. 


I ha , 2,ggd 

9o(— («xo> 


/ \ 

n X I -(ccrc) 


( 2 . 11 ) 


( 2 . 12 ) 


(we consider here only spatially homogeneous case, for simplicity), Aj° also satishes 


d,Af 


‘2gsd 

h 


n X Aj^. 


(2.13) 
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Substitution of eq. fl2.1ip into eq. fl2.13p leads to 



h + O((9o)^) = 0. 


( 2 . 14 ) 


Since this equation is true for an arbitrary n, we arrive at the relation = hSd 
Thus, we can obtain rris from Sc without any detailed calculation. The point is that the 
equation (12.131) is not the equation of motion of n, although it involves the time derivative 
of n. [The equation of motion of n is = 0, or, dQ{Sin + Aj^) = 0.] Before integrating 
out the electrons, it was the equation of motion of (ccrc) fl2.12p and, after integrating out 
the electrons, the equation fl2.13p determines the structure of AiSefr[n]. Conversely, the 
relation = hSc/i^g^s) must hold in order that the effective action AiSefj reproduces 
the equation of motion of {ccrc). We can repeat this procedure to arbitrary orders in the 
derivative expansion: hrst, write down all possible terms in the effective action AiSefr, with 
their coefficients left undetermined; second, derive the angular momentum Aj° from that 
AiSeflf and substitute it into eq. fl2.13p : then we can obtain the recursion relations between 
the coefficient^^ . For example, AiSeg to the fourth-order can be obtained as follows (see 
Appendix O for details): 


A5efr[n] =Sc 


d?x 


dt 


0(cos6> — 1) + 


1 
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2 

n ■ {h X n) + 


1 

16 



5 


-) 

9sdJ 


In 


2\2 


o(a; 


( 2 . 15 ) 


The same procedure can be applied also for the terms in the effective action AiSeg which 
involve the spatial derivative. 


3 Dynamical behavior of spin with inertia 

In this section, we describe classical dynamics of spin with inertia. For that purpose, it is 
helpful to use two equivalent pictures, which are summarized in Fig{T] One is a symmetric 

^Such recursion relations can be obtained also via the equation of motion of spin, as follows. The original 
equation of motion of localized spin is Sin = —SiB x n — JiSfdfn x n — gsd<A{c(Tc) x n. Substitution of 
the expression of (ccrc) [eg. (12.Ill) ] into this original equation of motion gives the term proportional to n, i.e. 
—gsdo.^ {ccrc) xn = — {2gsd'ms / fi)n+O {Oq) . The coefficient of this term is identical with the spin polarization 
Sc of the conduction electron, which gives the renormalization of the spin amplitude, Si —>■ Si+Sc- Therefore, 
we obtain the relation rus = hSc/{‘2gsd)- 
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Figure 1: The dynamics of spin with inertia, a symmetric spinning top, and a massive charged particle 
on a sphere subject to a monopole magnetic field Bm, are classically equivalent. They undergo precession 
motion accompanied by nutation under applied fields, which are a magnetic field B for spin, a gravitational 
field g for the top, and an electric field £ for the particle, respectively. See the main text for the detail. 


spinning top, and the other is a massive charged particle on a sphere subject to a monopole 
magnetic held. We consider here only spatially homogeneous spin, diti = 0, under time- 
independent magnetic held, for simplicity. 

3.1 Equivalence to a symmetric spinning top 

The equivalence between the classical dynamics of spin and a spinning top has been recog¬ 
nized in the literature, e.g. refs. na and [lO]. The content of this subsection is essentially 
a recapitulation of these facts, which we describe here in order for this paper to be self- 
contained. 

First, let us write down the Lagrangian of spin with hnite inertia and that of a spinning 
top. The Lagrangian of spin (12.Th is 

HI 

l^spin = + 5'0(cos 6^ — 1) -b BS cos 0, (3.1) 

while the Lagrangian of a spinning top in terms of the Euler angles (6*, 0, 'ip) is 

Ttop = 0COS dY -b figl cos9, (3.2) 

where Ji, Is are the principle moments of inertia, /i is the mass of the top, g is the gravita¬ 
tional acceleration constant, and I is the distance between the center of mass and the fixed 
extremity of the top. Here we have taken a symmetric spinning top and set two moments of 
inertia equal, Ji = l 2 - We take the positive directions of the external magnetic field B and 
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the gravity both in the positive z direction. Note that L^pin is a function of {9, 0) while L^op 
is that of Let us see below that they have equivalent dynamics concerning (6*,0). 

The equivalence can be directly seen at the level of their equations of motion. The 
equation of motion of spin is 


rris (O — 0^ sin 9 cos 9) + S^sm.9 + BSsm.9 


d 

dt 


m<i0sin^ 9 + S{cos9 — 1) 


0 , 


0 , 


while the equation of motion of the spinning top is 


(3,3) 


Ii{9 — 0^sin0cos0) + Is{ip + 0cos6*)0sin0 + figlsm9 = 0, 

d 


dt 

dt 


Ii4> sin^ 6 ^ + 13(0 + 0 cos 9) cos 9 


= 0 , 


-^ 3(0 + 0cos6*) 


= 0 . 


(3.4) 


From the last equation in eq. (13.411 . the canonical momentum M 3 conjugate to 0, 


M 3 = higjj + 0cos6'), 


(3.5) 


is conserved. Substituting this M 3 for 0 in the other two equations in (13.4p . we obtain the 
same equations as (13.3p with replacements 


Tfig i-y /i, S TFj, BS 


(3.6) 


Thus, the classical behaviors of 9 and 0 are the same for spin and a spinning top. 
We can see the correspondence more explicitly through their Hamiltonians : 

1 (M^ — S cos 9)'^ 


H Pg I 

■^spin — ^ “r 


and 


H - + 

tttop — ——r 


2m,, ' 2m, sin^ 9 

1 (M^-M 3 cos 0)2 


BS cos 9 


- fxglcos9 + 


(3.7) 


(3.8) 


2 Ji 2 Ji sin^ 9 

where pe = dL/d9 and M^ = dL/dcp are the canonical momenta of 9 and 0, respectively. 
These two Hamiltonians are completely the same under the replacements (13.6p . [The last 
term in (13.81) does not contribute to the dynamics of 9 and 0.] 

The Lagrangians (13.ip and (13.2p are related by the Legendre transformation about 0: 


Lspin(0,0,0 ; ^) = Ltop(0,0,0,0) - ^0 |p=p(0,0,s) (3-9) 

with replacements (13.6p . In the right hand side, 0 is substituted by S (or M 3 ) via (13.5p . 
The situation is quite similar to that of the familiar centrifugal force problem. There, the 


9 

























Figure 2: (a) Without any magnetic field, general motion of spin with finite inertia is a free precession 
motion, just like a spinning top. The spin Sn processes around the total angular momentum . (b) When 
a constant magnetic field B is applied, the total angular momentum processes around B. Therefore, the 
‘free precession cone’ in (a) processes around B as a whole, which corresponds to the Larmor precession in 
the absence of the inertia. What was called the free precession in (a) is now called the nutation. 


original Lagrangian is given as L{r,f,<J)) = (m/2)(r^ + — U{r), and we can obtain 

0-reduced Lagrangian by the Legendre transformation about <p\ M) = L — Mcj) = 

{m/2)P — M‘^/{2mr'^) — U{r). In exchange for reducing 0, there appears a hctitious potential 
M^/(2mr^). Likewise, spin dynamics is the '0-reduced dynamics of a spinning top, and the 
spin Berry phase (the second term in eq. fl3.ip ) appears as the hctitious potential arising 
from the '0-reduction. When we perform the Legendre transformation also about 6 and 0 
on both sides in eg. (13.91) . we are led to the same Hamiltonians (13.7p and (13.81) . 

Since the classical dynamics of spin with inertia and of a spinning top are equivalent, 
spin with inertia behaves in the same manner as a spinning top does (Figl2]). When a 
magnetic held is not applied, spin undergoes free precession: the spin precesses around the 
total angular momentum = Sn + rngn x h feq. (l2.10p i. which is a constant of motion. 
When a magnetic held is applied, spin undergoes the Larmor precession accompanied by the 
nutation: the total angular momentum precesses around the magnetic held (the Larmor 
precession), and the spin precesses around at each time (the nutation) [7j. The free 
precession solution of the equation of motion (13.3p (with .B = 0) is 

S 

0 = Oq (const.), 0 =-. (3.10) 

nis cos t^o 

Therefore, the free precession frequency, or the nutation frequency, cjq, is cjq ~ S/mg, 
assuming that cos 6*0 is order of unity (that is, the radius of the spin free precession or 
nutation is not very large). This gives the typical time scale for the inertial term, fmertia ~ 
mg/S. 
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Finally, we mention that usual spin with zero inertia, = 0, can be regarded as follows. 
It corresponds to the case of /i(= h) = 0 under replacements fl3.6p . This means that the 
other principle moment of inertia, I 3 , vanishes since, for an actual rigid body, any one of 
the principle moments of inertia is equal to or less than the sum of the other two, e.g., 
/3 < /l + /2 HE]. Thus, setting Ji = /2 = 0 leads to J 3 = 0. Therefore, there does not exist 
a spinning top corresponding to spin with zero inertia, in a non-relativistic framework. 


3.2 Equivalence to a massive charged particle 

Since the spin Berry phase corresponds to a monopole gauge held [1], the classical dynamics 
of spin is also equivalent to that of a charged particle on a sphere subject to a monopole 
background. Let us describe this equivalence and use it to understand the behavior of spin 
with inertia. 

A magnetic monopole yields a magnetic held in the radial direction with the strength 
inversely proportional to the square of the distance from the monopole. Being put at the 
origin x = 0 , the gauge held and the magnetic held of the monopole are 

. 1 — cos 6 * ^ „ 4 Q . N 

Q ; 7 ^ ^ (3.11) 

r sin d 

where x = r(sin 6 'cos(/), sin 0 sin(/), cos 6 *), e<^ is the unit vector in 0 -direction, and q is the 
magnetic charge of the monopole. The Lagrangian of a massive charged particle coupled to 
this monopole magnetic held and a constant electric held £ is 


L ~ $ 

TTi 

=—x‘^ + q^{cos 6 — l) + £-x, (3.12) 


where we take the electromagnetic scalar potential $ = —£ ■ x so that —V$ = £. This 
Lagrangian is identical with that of spin with inertia fl2.7p when we assume that the particle 
is constrained on a unit sphere and identify the direction of the particle x with the spin 
direction n, the mass m with the spin inertia the magnetic charge q with the spin 
amplitude S, and the electric held £ on the particle with the magnetic held SB on spin (we 
use calligraphic letters for electromagnetic helds on the particle). 

Being restricted on the unit sphere, the dynamical degrees of freedom of the particle are 
the spherical angles 6 and 0 , and the equation of motion is, in vectorial form. 


f d dL dL\ 

\dt 06 96 J 


f d dL dL\ 
sin 6 \dt 0 ^ dcj)) 


0 
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or mx = qxxx + £ — (£-x + miP‘)x^ (3.13) 

which is identical with fl2.9p (with diti = 0). The last term proportional to x in the right 
hand side of the second line is the constraint force, which keeps \x\ = 1. The particle is on a 
nnit sphere and is snbject to the monopole magnetic field Bm = qx (on the sphere \x\ = 1 ) 
and the electric field £. The general motion of the particle is, in the absence of the electric 
held £, the cyclotron motion, and, in the presence of the electric held £, the £ x Bm drift 
motioLL. 

In view of the eqnivalence between the dynamics of spin with inertia and of the charged 
particle on a sphere, the free precession of spin described in the last snbsection corresponds 
to the cyclotron motion of the particle with the freqnency Uq ^ q/m = S/mg, and the 
Larmor precession of spin accompanied by nntation corresponds to the £ x Bm drift motion. 


We snmmarize here the three pictnres for the classical dynamics described by (12.7p . 


pictnre 

inertia 

intrinsic magnetic held 

spin 

the inertia mg 

the amplitnde S 

charged particle 

the mass m 

the magnetic charge q 

symmetric spinning top 

the moment of inertia Ji 

the angnlar momentnm M 3 


We have seen in this section that spin with hnite inertia has an intrinsic precession mode 
(i.e. free precession or nntation). The magnetic held cansing this precession is the monopole 
magnetic held intrinsic to spin, i.e., the spin Berry cnrvatnre. The freqnency of the intrinsic 
precession mode, cjq ~ S/rris, is inhnite when the inertia of spin rUg is zero, bnt goes down 
to hnite valne when rrig becomes non-zero. 

As related works, we mention refs. IIZl[IHl[I 21 EniEIlE 2 ] where nntational motion of spin 
in Josephson jnnctions or tnnnel jnnctions between ferromagnets was discnssed by examining 
the Landan-Lifshitz-Gilbert eqnation with time-dependent Gilbert damping coefhcient. 

4 Behavior of spin with spatial inhomogeneity and un¬ 
der time-dependent magnetic field 

We have so far assnmed spatially homogeneons conhgnration of spin nnder time-independent 
magnetic held. Let ns inclnde spatial variation of spin and time-dependence of applied 

^Our convention for the sign of magnetic field is such that the equation of motion of a particle is mx = 
B X X + £, which yields drift motion of the particle in B x S direction, rather than £ x B direction. 
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magnetic field. The most typical behavior of spin under time-dependent magnetic held is 
the resonance phenomenon. In ref.[H], it is shown that spin with hnite inertia has a 
resonance peak at the intrinsic frequency cuq ~ S/mg. We discuss other behaviors of spin, 
such as spin wave, domain wall motion. We also discuss an unusual behavior of spin under 
a large and time-dependent magnetic held. 


4.1 Spin waves 

Let us consider spin wave solutions of spin with inertia, and see that there exists a new 
gapful spin wave mode with unusual handedness. When we assume n = {nx,ny,l) with 
Ux^riy <C 1, the equation of motion iQ from the action fl2.7p . 


m^n = Sn x n + JS^dfn 


(4.1) 


where we have set .6 = 0. Since the equation of motion is of second order in the time 
derivative, we have now two spin wave modes. When we substitute a plane wave -|- iriy = 
j^^-i{u}t-kx) equation of motion, with A an arbitrary small constant, we obtain the 

dispersion relation (FigJSj), 


2 rUs 


1 ± v^l -I- AJuiak?^ 


jse + 0{k^) 

-{i: + 


(4.2) 


One is the usual spin wave mode and the other is the new spin wave mode, the latter of 
which is essentially the free precession of spin mediated in space by the interaction J with 
neighboring spins. We have expanded the dispersion relation fl4.2p in powers of k and cut 
higher order terms because we are in the framework of the derivative expansion in deriving 

eq.([22D- 

The usual spin wave mode is clockwise seen from the positive ^(-direction, while the 
new spin wave mode is counterclockwise. In fact, to the leading order of the dispersion 
relation fl4.2p . the former obeys the equation of motion Sn = —Jdfn x = +JkSn x e^,, 
while the latter rngfi = —Sn x e^, where is the unit vector in ^-direction. Therefore, 
their handednesses are opposite. The total angular momentum in fl2.10p also rotates in 
counterclockwise for the new spin wave mode. 

When the spin and the angular momentum align in the same direction, = Sn, the fluc¬ 
tuation of spin around the stable conhguration is always clockwise. This can be understood 

^The assumption nx,ny <C 1 means that we consider the problem in the tangent plane on = 1, so the 
constraint force in (12.91) vanishes. 
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Figure 3: (a) The trajectory of a spin, which fluctuates around Sz = S point (this fluctuation propagates 
in space as spin wave). In the figure, two fluctuation modes are superimposed. One is the usual gapless 
mode rotating in clockwise seen from the positive z-direction (denoted by (2)), while the other is a new 
gapful mode rotating counterclockwise (denoted by Q). (b) The dispersion relation for the two modes of 
spin wave. Besides the usual gapless spin wave mode (in green color and denoted by (2)), spin with inertia 
has a gapful spin wave mode (in red color and denoted by ©), which is nutation of spin propagating in 
space by the spin-spin exchange. See eq. (IQ) . 


as follows (see e.g. ref. [23]). Let the stable configuration be = 1. The Poisson commuta¬ 
tion relation of the angular momenta, [j^,jy] = = S, leads to that of spin, [nx,ny] = 1/S. 

Therefore, the x and y components of spin are canonically conjugate to each other. More¬ 
over, the Hamiltonian of spin expanded around the stable point is H = + {ny)"^] 

with D > 0. These two facts, i.e. the canonical commutation relation of and Uy, and 
the form of the Hamiltonian, indicate that the fluctuation of spin around = 1 has the 
same dynamical structure as that of a one-dimensional harmonic oscillator with coordinate 
X and its canonical momentum p governed by the commutation relation [x^p] = 1 and the 
Hamiltonian H = D{x^ +P^)- The trajectory of the harmonic oscillator is always clockwise 
in the phase space. Therefore, the fluctuation of spin is also always clockwise. However, 
when the inertia of spin is hnite, the spin and the angular momentum generally do not align 
due to the additional contribution m^n x n in fl2.10l) . and moreover the Hamiltonian of 
spin such as eq. fl3.7p does not consist only of the potential term but also of the kinetic term. 
Therefore, the fluctuation of spin with inertia is not restricted to clockwise motion. 

4.2 Domain wall 

Let us see that a domain wall also has the free precession mode, corresponding to that of 
each spin. When easy-axis anisotropy (iLS'^/2) cos^ 6* is included to the action (12.7p . the 
system has a static domain wall solution, cos6 = tanh[( 2 ; — JT)/A], 0 = 0o (we take z-axis as 
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the wall normal) where X and 0o are integration constants, and A = JjK. Promoting X 
and 00 to dynamical variables, substituting the domain wall solution into the action 02.71) 
and integrating it over space, we obtain the Lagrangian of X and 00^ 

+ aV5) + - XM (4.3) 

with = rngN^/X"^ and = 2XA/a^ where A is the cross sectional area of the domain 
wall (we discarded irrelevant constant terms in the Lagrangian). We are going to show that 
this Lagrangian is of the same form as that of a charged particle on a cylinder with a constant 
magnetic held B = SN^/X"^ perpendicularly penetrating the surface of the cylinder. Take an 
orthogonal coordinate frame x = (x, y) on a cylinder, with y ~ y+2'KX the periodic direction. 
(The directions x, y and the direction outward normal to the surface of the cylinder make an 
orthogonal right-handed frame.) Then, the gauge potential A. for the constant magnetic held 
B in the outward normal direction of the surface can be taken as Ax = By/2, Ay = —Bx/2, 
and the Lagrangian of the particle is 

TD 7H Lj 

L = -P-A-x = -(x^ + y^) + -{xy- xy) (4.4) 

Therefore, when we identify (X, A0o) with {x,y), the Lagrangian fl4.3p is that of a charged 
particle under a perpendicular magnetic held of hux B = SN^/X"^ (FigJl]). The magnetic 
held B for the particle originates in the spin Berry phase, and the mass of the particle 
comes from the inertia of spin in the action fl2.7p . 

As a result, the particle undergoes cyclotron motion (in other words, the domain wall un¬ 
dergoes the intrinsic free precession mode), where X and A0o oscillate, X = r cos{uot), A0o = 
rsin(a;ot) with r a constant, at frequency uq = B/M„ = SN^/{X‘^M^) = S/rris- This fre¬ 
quency corresponds to that of spin free precession discussed in section [3l When an external 
held is applied on this domain wall, the particle on a cylinder {X, A0o) feels the correspond¬ 
ing force on it, and it moves in the direction perpendicular to the force, accompanied by the 
cyclotron motion (nutation). That is, the particle undergoes £ x J3 drift motion, where £ 
is the force and B is the magnetic held penetrating the cylinder {X, A0o) perpendicularly. 

The discussion above is well valid when the inertia rris is so large that the intrinsic 
frequency Uq ^ S/rris is less than the frequency of spin waves: S/rUs < KS'^/h. When S/nis 
exceeds this frequency scale, the dynamics of the domain wall can be no longer described only 
by the collective coordinates X and 0o, and spin wave excitations should also be included 
into the dynamics. 
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A2 


Figure 4: The collective coordinates {X, X(j)o) of a domain wall can be regarded as coordinates of a particle 
on a cylinder. This particle is subject to a magnetic field fiux B = SN^^/X'^ penetrating the cylinder 
perpendicularly. The inertia of spin introduces the mass = nisN^/X'^ into this particle. 


4.3 Parallel shift of spin under time-dependent magnetic field 


In the spin nutation such as in FiglU spin moves up and down parallelly to the given magnetic 
held, but this motion is oscillatory and the guiding center (i.e. the center of the nutation) 
moves only perpendicularly to the magnetic held, just as the usual Larmor precession. Is it 
possible that the guiding center motion becomes diherent from the usual Larmor precession 
and has a hnite component of velocity parallel to the magnetic held, considering that the 
spin with inertia does not obey the usual Landau-Lifshitz equation? In this subsection, we 
see that when the magnetic held is time-dependent, the guiding center motion has hnite 
component of velocity parallel to the time-derivative of the magnetic held. 

Bearing the particle analogy discussed in section [3] in mind, this parallel motion is most 
directly seen by looking at the behavior of a charged particle, located on a plane (x, y) under 
a constant perpendicular magnetic held B and a time-dependent in-plane electric held S {t) 
(recall that an electric held S on the particle corresponds to a magnetic held B on spin). 
The general solutions of the equation of motion is given by {x,y) = (^, 77 ) -b {X,Y), such 


that 


and 



(4.5) 



B1 + u-‘‘{<Pldfi) 


u^^d/dt —1 
1 uj~^d/dt 



B 

OJn = -. 

m 


(4.6) 


That is, we can divide the motion into the cyclotron motion (,^, rj) and the guiding center 


motion (X, Y) even when the electric held is time- 


!-dependeni@. 


We can see that the guiding 


^As long as ^d/dt < 1 and the operator (1 -|- ‘^d^/dB) ^ is well-defined. 
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center motion {X, Y) has a finite component of velocity parallel to the electric field, when the 
mass is finite and the electric field is time-dependent. As we saw in section [31 the classical 
dynamics of spin with inertia is a spherical version of this planar problem, with the role of an 
electric field for the particle played by a magnetic field for spin, and therefore the essential 
features will be the same (especially when the nutation radius is small and the curvature 
of the sphere can be neglected): (i) the spin nutation and the guiding center motion are 
decoupled; (ii) the effect of the inertia on the guiding center motion is significant when the 
time rate of change of the applied field is comparable to the intrinsic frequency S/mg. 

Therefore, the effect of the inertia is not only the superimposition of the nutational 
motion to the Larmor precession. The inertia of spin gives an independent effect on the 
trajectory of the Larmor precession itself, cooperatively with the time-dependence of the 
applied field. Even when the nutation radius is small and there are no signs of nutation in 
experimental data about the time-profile of magnetization, it never indicates that the inertia 
has no effects on the dynamics, especially when the applied magnetic field has very rapid 
time-dependence, such as ultrafast magnetization. 

We illustrate this fact by numerical calculation (FigJS]). We fix a magnetic field in the 
2 ;-direction. First, we make the magnetic field in the negative 2 ;-direction, then change it 
into the positive ^-direction. We can see that the spin moves parallelly (i.e. in the positive 
^-direction) to the time derivative of the magnetic field when the magnetic field is changing. 
This phenomenon is similar to usual spin flip by the Gilbert damping, but is different in that 
we do not have to vary the direction of the magnetic field and moreover the time change 
rate of the magnetic field can be arbitrary (as long as our derivative expansion is valid) and 
is irrelevant to the damping coefficient. 

This unusual behavior of spin can be simply understood when we look at the angular 
momentum rather than spin itself. As in fl2.10p . the angular momentum of spin with finite 
inertia is no longer Sn but has additional non-adiabatic contribution m^n x n. When 
the Larmor frequency B is much smaller than the intrinsic frequency (nutation frequency), 
B S/rUg, then the vector m^n x n points inward to the center of the nutation circle, and 
the trajectory of the angular momentum of spin is overlapped with that of spin (Fig[ni(a)). 
On the other hand, when B ~ S/mg, nutation does not form circles but rather ripples on 
the trajectory of spin. In this case, the vector x n points in almost the same direction 
at each time. As the result, the trajectory of the angular momentum deviates from that 

®The switching of the sign of the magnetic field from negative to positive is not essential for the shift of 
the spin in z-direction. Only the time-dependence of the magnetic field is essential. 
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Figure 5: Response of spin with finite inertia to a time-dependent magnetic field B{t) = Bq tanh(/(t — to)) 
in z-direction with Bq = lOOTHz and / = 27rTHz = 6.3THz. We take the spin amplitude a.s S = h and the 
inertia as rris/h = 0.1ps/(27r) = O.Olfips. (a) Time profile of Sz/S{= riz) of spin and applied magnetic field 
B. (b) The trajectory of spin. The time is divided into three domains I-III. In I, the magnetic field points 
in the negative z-direction and spin processes counterclockwise seen from the positive z-direction. In II, the 
magnetic field switches to positive z-direction and the z-component of spin, Sz/S, also changes accordingly. 
In III, the magnetic field points in positive z-direction and spin processes clockwise. 
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of spin (FigEl^b)). FigJT] is a replot of the trajectory of spin in FigJS] together with that 
of the angular momentum. The vector m^n x n points in the positive 2 ;-direction in the 
time domain I (described in FigjS]), and in the negative 2 ;-direction in the time domain III, 
resulting in the conservation of the ^-component of the angular momentum. 

This parallel shift of spin in the direction of the time derivative of the magnetic held can 
be understood also from a viewpoint of a spinning top, where the magnetic held on spin 
corresponds to the gravitational held on the top. The time-variation of the gravitational held 
can be realized hctitiously by putting the top in an elevator accelerated in the negative 2 ;- 
direction. There, the top feels, besides the real gravitational held in the negative ^-direction, 
a hctitious inertial force in the positive ^-direction due to the acceleration. This inertial force 
in the positive 2 ;-direction makes the top stand up in the positive 2 ;-direction, which explains 
the shift of spin in ^-direction. Thus, the inertia of spin yields an inertial force, literally. 
In this respect, the connection between the inertia of spin and the inertial ehect on spin 
[2n [25] will also be interesting. 

5 Summary and discussion 

We examined the inertia of spin in metallic ferromagnets, which arises in the derivative 
expansion of spin ehective action and represents non-adiabatic contribution from environ¬ 
mental degrees of freedom. We derived an explicit expression of spin inertia in an sd model. 
The equivalence between the dynamics of spin with inertia, spinning top and a charged 
particle on a sphere was explained. The behavior of spin with spatially inhomogeneous con- 
hguration or under time-dependent magnetic held is studied. The hnite inertia has mainly 
two effects: the superimposition of the nutation on the usual Larmor precession and the 
parallel shift of the Larmor precession. 

As explained several times in this paper, most of unusual behavior of spin with hnite 
inertia originates in the non-adiabatic component m^n x n in the relation between spin 
and the angular momentum, = Sn -|- nisTi x n. Due to this diherence between them, 
spin acquires a kind of separate freedom from the angular momentum, whose behavior is 
relatively restricted due to the conservation law. As the result, spin undergoes free precession 
or nutation discussed in section [3l and furthermore parallel shift along the direction of the 
time derivative of applied magnetic helds, discussed in section |U Since one of the main 
purposes of spintronics is understanding and application of diverse cooperative phenomena 
between localized spins and their environmental degrees of freedom, researches beyond the 
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(a) i? <C S/llls 



trajectory of 

-spin trajectory angular momentum 

-- m.sn X h 

Figure 6: Trajectories of spin Sn and angular mo¬ 
mentum = Sn+irisnxri. (a) When B <C S/nis, 
spin nutation forms circle shape. The additional 
component rngn x n points inward to the center 
of the nutation, and the trajectory of the angu¬ 
lar momentum sits within that of spin, (b) When 
B ^ S/tUs, spin nutation forms ripple shape. The 
additional component nisTi x h points to one side 
of spin trajectory. As the result, the trajectories of 
spin and the angular momentum deviate. 



Figure 7 : Replot of the trajectory of spin in FiglS] 
(from slightly a different view angle), together with 
that of the angular momentum. The magnetic field 
is R ~ lOOTHz and S/m^ - 207rTHz = 62.8THz, 
corresponding to the case (b) in Figl6l Although 
the z-component of spin changes in time, that of 
the angular momentum is constant in motion. The 
axes are in unit of S. 
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adiabatic limit will be well-motivated for this spirit of spintronics; in the adiabatic limit, the 
environmental degrees of freedom are rather obedient to the localized spins. In this paper, we 
estimated that in metallic ferromagnets such non-adiabaticity is effective for sub-picosecond 
dynamics, i.e. ultrafast magnetization. For broader applications, physical situations where 
non-adiabaticity is effective for sub-nanosecond dynamics are also desirable. 

Apart from the context of effective dynamics (i.e. the dynamics including the effects of 
environmental degrees of freedom), the inertial term, together with the spin Berry phase 
term, can be present in the dynamics of the order parameters in systems whose ground 
states can be intermediate states between ferromagnet and antiferromagnet (see e.g. |26jb 
The discussion of this paper will be also applicable to such systems. 
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A spin effective action 

In this appendix, we calculate the spin effective action by integrating out a conduction 
electron coupled to spin through the sd interaction, and derive the expression of spin inertia 
fl2.4p . We always make summation over repeated indices unless otherwise stated. 

The action of a conduction electron, represented by held operators c and c, coupled to 
localized spin n is 




with gsd the sd coupling constant. What we calculate in this appendix is the contribution 
to spin effective action from the conduction electron, AiSeff[n], dehned as 



(A.2) 


Adding this AiSeff[n] to the original spin action 



(A.3) 
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where 6 and 0 are the polar angles of the nnit vector n, we obtain the total spin effective 
action 5efr [n] = Ss [n] + AiSgff [n]. 

First, we diagonalize the sd coupling in flA.lj) . In order for that purpose, it is conve¬ 
nient and essential to employ the unitary-transformed electron held a such that c{x, t) = 
U{x,t)a{x,t). Here, the SU{2) matrix U is dehned so that 


U\n ■ cr)U = (Ts. 


(A.4) 


As such, U has indehniteness about the right U{1) gauge transformation, U ^ U exp(zx(a 3 , t)a‘i) 
with an arbitrary function. [Due to the dehning relation c = Ua of a, this unitary 

transformation corresponds to that of a, i.e. a —)■ exp{—ix{x,t)a 3 )a. The original electron 
held c remains invariant.] One choice of t/ is 17 = m cr with m = (sin | cos 0, sin | sin 0, cos |). 
Geometrically, U rotates the electron spin in the positive ^-direction, | '(), to the direction 
of n. For the choice U = m cr above, U\ t) = I^t-) with |n) = cos(6'/2)| 'I') -f-e*'^sin(6*/2)| 0). 
The indehniteness of U, U U exp(ix(T 3 ), corresponds to the phase rotation of |n), i.e. 

|n) —)■ exp(ix)|n). The sd coupling becomes diagonal in the a-frame, c(n • cr)c = da^a but, 
in compensation for the SU{2) gauge transformation by U, there appears a spin gauge held 
Afj, = —iWdfj,U = (m X ■ cr in the derivative terms, d^c = U{dfj, + iA^)a. 

Then, the action flA.lD becomes [27] 


Se= d'^xdt a I ih{dt + iAt) + 


lT?{di + iAi 
2 m 


+ ep + S'sdO's ® 


=iSo[a, a] -I- / d^xdt A^[—haa°‘a) + A“ ( ^—aad^ia ] -1- A“A“ 



(A.5) 


Here Sq [a, a] is 




5o[a, a] = / d'^xdt a ( ihdt + ^ ^ + ep + S'sd^'a ) a. 


(A.6) 


and A^ = A“cr“ (a = 1,2,3) is the spin gauge held (make distinction between the electron 
held a and the SU{2) index “), which is, in terms of the polar angles of n, 


A, 


^—dfj^O sin 0 — sin 6 cos (t)d^(j^ 
0^6 cos 0 — sin 6 sin 00^0 
y (l-cos6')0^0 y 


(A.7) 


where (A^)“ = A“. Because A^ is proportional to the hrst order derivative of n, the 
derivative expansion of AiSes in powers of corresponds to that of A^\ 


i 

h 


A5efr[n] 


i 


h 


d'^x 


A“(-h(aa“a)) + A“ 
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2\h 


d^xd^x' 


A1{x)A\{x')h^{aa°'a{x)aa^a{x')) 


3 ' 


+ A1{x)A^-{x') ( -—j {da^ di a{x)da^ dj a{x')) 


+ A'^{x)A'l{x') ( (acr“a(x)acr^^a(a;')) + 0{AA 

\ m J J 

where time-ordered quantum expectation values are taken about So, 

(...)= [ VdVa (...)eH‘5oM]^ 


(A.8) 


(A.9) 


and we take only the connected diagrams to evaluate the quantum expectation values. Be- 

^ _ y ^ _ y 

cause the non-perturbed action Sq flA.bp is isotropic in space, {da°‘ d io) and {da°‘a{x)da^ di a{x')) 


vanish. We dehne the notation for correlation functions as 






p°'°{x,x') = i{aa°'a{x)aa°a{x')), = i{aa^ di a{x)aa° dj a{x')). 


(A.IO) 


Using the Fourier components, we expand these non-local quantities. For example. 


p‘^\x,x') = 


d^q 


^iq-(x—x') ab 


p^\q) 


1 

=5^^\x - x')p“\=o A - ^')^l.=o + ■ ■ ■ . (A.ll) 


When we substitute this into the expansion flA.Sp . the terms other than the hrst term in 
flA.lip increase the order of the derivative on n after partial integration, e.g. A^{x)A\{x')d^5‘^^^ {x— 
x') —^ —d^[A^{x)A\{x')]5^‘^\x — x'). Therefore, we take only the hrst term in flA.lip as long 
as we are interested in at most (9((9^n)^-terms. After all, AiSefr[n] is, to (9^n)^-order, 

A5.*[n) = j d‘x (-h(an‘‘a))A1 + A|l{aa)i.,iy - ^Re4|,.d A^A] + AtA\ 

(A.12) 

(the imaginary parts of yA and contribute to damping and are neglected here). We 
calculate each term in the following. 

Firstly, we calculate (acr“a). The time-ordered Green’s function for So flA.6D is 

d^kdu 


-i{aa{x)ap{x')) =5ap 


(2vr) 


_^ik {x x') iui(t (no summation over a). 


9k,uj,a — I • r ’ 

CJ ^koi A 

(the indices a, (3, ■ ■ ■ = ± correspond to spin up or down) with 


hxxk± = 


2 m 


tF T a,d, V = sgn(uita) X 0, 


(A.13) 


(A.14) 
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where 0 denotes a positive infinitesimal. We define a 2 x 2 matrix by {gk,uj)ai 3 = ^ai 39 k,uj,a 
(no summation over a), or, 

9k,UJ = '^{.9k,txj,+ 9k,LU,— ^ {.9k,LU,+ 9k,u),— ))- (A.15) 


Then, 


(aa^a) = -i j = -iS^ X^(±) 


cPkdu 


^iujO 


(27r) 

=4“E(±) / 


(27r)'^ u — ujk± + ih 


± 


di^k 


9{-Uk±) = 


I^F+ ^F- ra ^F± _ I 

■ 03 , ^_ — cp ± g^A) 




2m 


(A.16) 


where 9{x) = 1 for a: > 0 and zero otherwise. 

Secondly, the number density of the conducting electrons n = {da) is given by n = 
(fc|+ +A;|._)/(67r2). 

Thirdly, we calculate ReXij\q=o- Using 

tr[(T 9k-\~q,,LO-^Q] ^ ^ [(^afo Sa3^b3 T '^^ab3^9k,LU,:k9k~^q,uj-\-Q,^ T ^a3^b39k,(jj,:L9k~\~q,LO-^Q,,^^ 

(A.17) 

and 


• f duj _ d(uJkoi)d{^ ^k-^qp) 

^ / ~^^9k^L0,a9k-\-q,0J-\-O,,l3 


^ka)d{^^k~\-ql3^ 


^k~\-q^ ^ka ^ T ^k-[-ql3 ^ka ^ 


(A.18) 


we obtain 

X?‘(q,n)=-i 

= E 


hjduj 

- :^{2ki + qi){2kj + qj)tT[a‘^gk,ojCT^gk+q,ui+n] 


(2f 

d?k 


3 {2ki + qi){2kj + qj) 


(27r 

X \ {dab — da3db3 T 
+ da3db3 


d{uJk:k)d{ CJ/g|g^) 


d{ ^k:k')d{(-^k~^q^') 


^fcih ^ d” ^d CU/g-|_qfzp a^fcdi ^ ^d_ 

d{uJk^^d{ UJk-^qzt) d{ (■^kzL)d {uJk-\-q^') 


_^fcH-qdi ^fcdi ^ T ^d a^fc+qdi ^fcdi ^ ^d 


(A.19) 


We are interested here only in the real part, Rex“-. Using 


—^— = P—=Fi7 r(5 (x) and 9{x)9{—y)—9{—x)9{y) = 9{x) — 9{y) = 9{—y)—9{—x), (A.20) 

x ± iO X 

then. 


Rex°^(q,n) =J2 / i0;(2*:. + 9i)(2% + fc) 
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.. J/'X XX I X X -D^(“^fc+9±) - ^(-l^fc±) 

X (daft - dasdftsjP—---- 7 ;;-h dasdbsP- 


^fc+g=F ^fc± ^ 


^k+qi: ^fc± ^ 


± £ab3Tr[0i-UJk+qT) - 0{-UJk±)]5{Uk+q^ “ UJk± “ 


(A.21) 


Let us calculate the value of each component of Rex“j'|g=o- First, 


r a3u k 

Rex]j\q=o = Rex^if\q=o =J2 J 7^4/ci/cj(±);^[0(-a;fc^) - e{-Uk±)] 


( 2 vr)= 

h 2(fc| 


'F+ 


U5 ' 
Fp-; 


9sd 


IStt^ 


25 'sd 


■6ij. 


(A.22) 


These components do not depend on the order of taking the g —)■ 0 and hi —)■ 0 limit. On 
the other hand, the value of the component Rexff depends on the order: From eq. flA.2ip . 
it is clear that when we take the q —)■ 0 limit hrst, Rexfj vanishes, while when we take the 
O —)■ 0 limit hrst and then take the q —)■ 0 limit, it gives 0/0 and, using the L’Hopital’s rule. 


lim lim Revf/ 
q^o n^o ^ 




4mn 


(A.23) 


The order of the q ^ 0, O —)■ 0 limits can be simply determined by the isotropy of the 
internal space (i.e. the space where the vector n is dehned). For example, the combination 
+ {AjY = (l/4)(c}jn)^ (see eq. flA.7p i is isotropic (i.e. does not depend on the rotation 
n Rn with a constant SO{3) matrix R), but {A^Y = (1/4)(1 — cos0)^(c?j0)^ is not an 
isotropic quantity. When we employ hmn_,.o limq_j.o Reyf^ = 0 as the value of Rexf||g =05 
then the (Af)^-term appears in AiSefr flA.12p due to the nonzero (aa)-term. On the other 
hand, when we employ eq. flA.23p as the value of Rexf||g=o, then the (A^)^-term vanishes i 


m 


A5efr due to the cancehation with the (aa)-term. Therefore, we should employ eq. flA.23p as 
the value of Rex/j |g=o 




For non-diagonal elements, it is clear that ReXij\q=o = 0 regardless of the order of the 
q —)■ 0 and O —)■ 0 limits. 


^ The order of q —>■ —>■ 0 limits can be determined also by the t/(l) gauge symmetry U —>• 

U exp(ixCT 3 ), a —>■ exp(— 1 x 173)0 described below ea. (IA.4l) . Under this gauge transformation, the spin gauge 
field Afj^ = A“cr“ changes as 


^ /^cos2x -sin2x\ /AA 
\Al) Vsin2x cos2x ) \Al) ’ 


Ajl —>■ + d^x- 


(A.24) 


The effective action A5eff (1A.12I) should be invariant under this gauge transformation. Therefore, (Aj)'^ + 
(Af)^-term is allowed but (A^)^-term is not. 
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Finally, let us calculate Rep“^|g=o- The calculation is almost the same as that of Rex“j'|g=o: 
just drop the factor —{2ki + qi){2kj + qj) from eq. flA.19p . Then, 


Rep^\=o = Rep^\=o = -J2j “ ^(-^fc±)] 

h k^_^_ k^_ 

gsd 67r2 


(A.25) 


The order of taking the q, —?■ 0 limits matters again for Rep^^. This time, hmq_j.o hmQ_j.o Rep^^ 7^ 

0 while limf7_j.o limq_,.o Rep^^ = 0. Again, we determine the order of the limits by the 
isotropy about n. When Rep^^|g=o does not vanish, appears in AiSes, which breaks the 
isotropy. Therefore, contrary to the case of Rex^^^|g=0) we employ limf7_j.o hmq_j.o Rep^^ = 0 
as Rep^^|q=o- Finally, Rep^^|q=o obviously vanishes. 

After all, the expression of the contribution from the conduction electron to spin effective 
action is 


AiSpff = / d^x 


d?x 


p3 

- 




kl_ f 1 k% 

At + 


F+ 


p5 

f^JT— 


2m \ 2m Agd ISvr^ 


- n) 1(4')= + (4)=1 


fc2 p3 _ p3 

^-[{Air+iAtr] 


2A 


sd 


bvr^ 


dt 


Sc^{cose - 1) - 


(A.26) 


where a is the lattice constant and we have dehned the spin polarization Sc of the conduction 
electron per lattice site, the spin-spin exchange coupling Jc and the inertia of spin as 


e_3 ^ ^'f+ ~ ^F- 

2 67r2 


aS , 


fe2 1 p5 _ 

It 1 

m Qsd 307r2 


m., = 


hSc 

2psd’ 


respectively. It is easy to see that Jc > 0 by rewriting it as 


Jr = 


{kp^ — kpJf 


(fcF+ + ^kp+kp^ + kp_). 


(A.27) 


(A.28) 


5'2 1207r2m kp^ + kp_ 

We may safely neglect the effect of impurity, because the lifetime r of the conduction elec¬ 
trons will enter only as h/{epT) or h/(psd'r), both of which we assume much smaller than 
one. 

Adding this AiSeff[n] to the action of localized spin fl2.6p . we obtain the total effective 
action iSefr[n], 


Sesln] = 


d?x 


dt 


S'0(cos 9 — 1) 


/n \2 I '^s . 2 

-^(diu) +^n 


(A.29) 


with S = Sc + Si and J= {JcSl + JiSf)/S^ 
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B Derivation of the angular momentum (l2.1()l ) 


In this appendix, we describe the derivation of the angular momentum current j*) fl2.10l) 
from the effective action iSefr fl2.7p . for the readers who ar^not familiar with such derivation. 

lis rotation, a vector v 


Consider a SO(3) rotation in internal space (spin spacejj. Under t 


I 


When the rotation is 


in internal space, such as spin n and magnetic held B, is rotatec 
inhnitesimally small, it is given by 6v°'{x,t) = where 0“ are constant inhnites- 

imal transformation parameters. The effective action (12.7p is invariant under this rotation. 
The Noether current corresponding to this invariance, i.e. the angular momentum current, 
can be obtained by the following trick. Make the transformation parameters be depen¬ 
dent on the position in the real space and time, 6°‘{x,t). Then, the variation of the action 
under this space-time dependent internal rotation does not vanish, but is given in the form. 




d^x 


dt{doO ■ f + diO ■ f) = - 


d?x 


dt 6 ■ {dof + dif). 


(B.l) 


Here, terms without any time- and spatial-derivatives on are absent, because the rotation 
would be a symmetry of the action if 6*“ were constants. We have performed the partial 
integration in the last equality in eq. flB.H . When n satishes its equation of motion, then 
hiSeff = 0 because such n gives a saddle point of iSeg. Since hiSeg = 0 for arbitrarily position- 
dependent 6°‘{x,t), it follows that + dij^ = 0 when n satishes its equation of motion. 
Thus, (j°, J*) is the conserved current corresponding to the symmetry under internal rota¬ 
tions, i.e. the angular momentum current. [In fact, since n is a unit vector, an arbitrary 
variation of n is given by a rotation in internal space. Therefore, from the variation fIB.ip . 
the conservation law doj^ + = 0 is the equation of motion itself.] 

Concretely, the variation of each term in the effective action fl2.7p is given as follows. For 
the spin Berry phase term. 


d^x 


d?x 


d / ——dt S(j){cos6 — 1) =(5 / ——dt 


= [ ^^-^dt [ du S[du{6 ■ dtn 




du Sh ■ {dufi X dtfi) 

dt{6 ■ duu)] 


'0 


®An element R of 50(3) group and an element U of SU{2) group are related by Rab<S’ = with 

cr“ the Pauli matrices. We are interested in the classical dynamics of a spin vector, that is, the behavior of 
n = (^jcrl'b), where |^) is a quantum spin state. A 517(2) rotation on IiIj) causes a 50(3) rotation on n, 
i.e. {'4’\Wa°'U\'4’) = Rab'nt’- 

®The magnetic field B forms a scalar product with n in the Zeeman term, which means that B is also a 
vector in internal space with the same transformation property as n, within the framework presented here. 
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(B.2) 


= [fEjte-Sn, 

J O'* 

where we have temporarily used n{x, t, u) which is dehned over extended space-time (a:;, t, u) 
with 0 < M < 1 a dummy direction and with boundary conditions n(a;, t, n = 0) = n{x, t) 
and n{x,t,u = 1) = const. The transformation parameters 6 are also extended to 6{x,t,u) 
similarly. In the last equality in eq. flB.2p . we have performed the partial integration with 
respect to t. For the inertial term and the exchange coupling term, 


5 


5 

d^x 


d?x , rus . n 
—dt —-n? = 
a3 2 


d?x , X , . s 

——dt 0 ■ (nisn X n), 


dt 


-JS^ 


{diTif = 


d^x 


dt diO ■ {JS’^di'n x n). 


(B.3) 


The Zeeman term is irrelevant for deriving the angular momentum because it does not 
involve the derivative on n, which means that it does not yield terms proportional to d^O 
or di6 in eq. fIB.ip under the rotational transformation. In all, from eq. llB.ip . the angular 
momentum current can be read as 


= Sn -|- m^n x n, = JS’^diU x n. 


(B.4) 


C Derivation of higher order terms in spin effective 
action 


In this appendix, we derive higher order terms in spin effective action AiSefr (12.1 5p in the 
way discussed in the last part of section |2l We concentrate only on time-derivative terms, 
but derivation of higher order terms with spatial-derivatives is essentially the same. First, 
all possible terms in spin effective action to the fourth-order are 


AiSpflF — 


-dt 




Dc0(cos 6^ — 1) 


(C.l) 


The angular momentum Aj° can be derived from this action in the way described in Ap¬ 
pendix [Bl as 


Af = ScTi m^n X n -|- C3(2n -|- 3n^n) -|- 2c4(n x fi — n x n) -|- Ac^n^n x n. (C.2) 


This must satisfy doAj^ {2gsd/h)n x Aj° = 0. Substitution of eq. flC.2D makes, with 

^sd = Qsd/ h, 

doAf + {2g^d/K)n x Af 
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{Sc - 2 usdms)n + (m* + 4wsdC3)n x n 

+ (2c 3 + 4a;sdC4)n + (6C3 + 12a;sdC4)(n • n)n + (803 — 4a;sdC4 — 8a;sdC4)n^n. (C. 3 ) 


The unique solution in order that each coefficient vanishes is that 

Sc rus C 3 „ 5c3 

'^s 7 ^ ) C3 - , C4 — , C4 — . 

ZCUsd ^l^sd O^sd 


(C.4) 
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